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(57) Abstract: The invention concerns the 
modelling of interactions between an incident 
wave and an obstacle, in particular in the field 
of non-destructive control. The invention is 
characterized in that it consists in: meshing 
the obstacle surface and assigning at least 
one source (S { ) to each surface element 
(dSj); then calculating conditions at the 
boundaries in each mesh of the obstacle and 
deducing therefrom source values; based on 
an interaction matrix and said source values, 
estimating a physical quantity representing the 
interaction between the wave and the obstacle, 
in any point of space. 

(57) Abrege : L invention concerne la modeli- 
sation des interactions entre une onde incidente 
et un obstacle, notamment dans le domaine du 
contrdle non destrucuf. Selon rinvention, on 
maille la surface de Tobstacle et Ton attribue 
au moins une source (Si) a chaque element de surface (dSi). On calcule ensuite des conditions aux limites en chaque maille de 
V obstacle et Ton en d6duit des valeurs de source. A partir d*une matrice d 'interaction et de ces valeurs de sources, on estime une 
grandeur physique representative de V interaction entre i'onde et Tobstacle, en tout point de Tespace. 
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Precede pour evaluer une grandeur physique representative 
d'une interaction entre une onde et un obstacle 

L 1 invention concerne la modelisation des interactions 
5 entre une onde incident e et un obstacle de cette onde, 
notatranent dans le domaine du contrdle non destructif. 

La modelisation des interactions entre une onde et un 
obstacle recevant cette onde, tel qu'une cible plac^e dans 
10 la zone sensible d'un capteur, trouve une application 
avantageuse en contr61e non destructif . 

On connait un proc€d§ de modelisation dit "par Elements 
finis" consistant & appliquer un pavage de l'espace 
15 tridimensional entourant 1' obstacle et a §valuer les 
interactions pr€cit€es pour tous les pav£s de l'espace. 

Les procedSs de calcul par "elements finis" apportent une 
resolution a un probleme pos§ sous forme d' Equations 

20 dif f Srentielles partielles. lis se fondent sur une 
representation de l'espace d'€tude par un assemblage 
d' elements finis, & 1'interieur desquels des fonctions 
d' approximations sont d^terminees en termes de valeurs 
nodales de la grandeur physique recherchee. Le problSme 

25 physique continu devient done un probleme discret aux 
616ments finis oil les grandeurs nodales sont les nouvelles 
inconnues. De tels proc6d6s cherchent done ai approcher la 
solution globale, plutSt que les equations aux d€riv6es 
spatiales partielles de depart. 



30 
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La discretisation de l'espace pris en compte assure que ce 
dernier soit entierement recouvert par des elements finis 
(lignes, surfaces ou volumes), cette operation porte le 
nom de "maillage" dans l'espace bidimensionnel (2D) ou de 
5 "pavage" dans l'espace tridimensionnel (3D). Les elements 
mis en jeu sont soit rectangulaires ou triangulaires en 
2D, soit parallel6pip6diques ou t6tra#driques en 3D, lis 
peuvent §tre de tailles differentes, repartis uniform£ment 
sur la surface ou non. 

10 

En general, la grandeur physique recherch€e, telle qu'un 
potentiel electrostatique ou une valeur de pression, est 
connue sur la frontiere du domaine. Cette frontiSre peut 
§tre fictive. On y impose des conditions aux limites. Le 
15 potentiel est done inconnu i. l'interieur du meme domaine. 
On d£finit alors un nceud comme etant un sommet d'un 
element. Les inconnues du probldme sont done les valeurs 
du potentiel en chaque noeud de 1' ensemble du domaine. 

20 A titre d' illustration, la figure 6 de l«art anterieur 
reprgsente un exemple de surface, constitute par deux 
materiaux Ml et M2, de proprigtes electromagnet iques 
differentes, et maillee par des elements triangulaires 
comportant chacun trois nceuds Ai, Bi et Ci . L ' ensemble du 

25 domaine est d€limite par une frontiere F. 

Une fois le maillage d^fini, plusieurs approches existent 
pour transformer la formulation physique du probl^me en 
une modelisation discrete aux elements finis. Si le 
30 probleme est formula par des Equations dif f 6rentielles et 
consiste & minimiser une f onctionnelle, on applique 
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genSralement une methode variationnelle . Cette 
transformation aboutit a une formulation matricielle dont 
la resolution donne les solutions nodales, les solutions 
aux points non nodaux etant obtenues par interpolation 
5 lineaire. 

Neanmoins, de tels calculs, en trois dimensions, 
reguigrent des ressources inf ormatiques importantes et 
ggnerent des temps de calcul tres longs, en depit de 
10 l'accroissement des performances de logiciels permettant 
d' impiementer ces calculs. 

Certes, les problemes 2D, souvent simplifies par des 
conditions de symetrie avantageuses pour ne modeiiser 

15 qu'une partie de la geom^trie, se resolvent rapidement. 
Mais il n'en est rien pour les problemes 3D, qui sont les 
plus frequents. La figure 6 montre combien la finesse du 
maillage, c'est a dire le rapport entre la taille d'un 
Element et celle du plus petit detail du domaine, fait 

20 croitre le nombre de noeuds. 

En consequence, la quantity d' equations et d # inconnues 
augmente proportionnellement , et, de 1&, le temps de 
calcul ngcessaire §l la resolution du probleme. Il est 
25 important de prgciser que la generation du maillage, a 
savoir la discretisation de l'espace de travail, et la 
generation de la liste des noeuds consomme un temps de 
calcul superieur & celui necessaire a la resolution du 
probleme. 

30 



La presente invention vient ameiiorer la situation. 
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Elle propose a cet effet un procede pour evaluer une 
grandeur physique associee a une interaction entre une 
onde et uri obstacle, dans une region de 1' espace 
tridimensionnel, dans lequel : 

a) on determine, par maillage, une pluralite 
d'echantillons de surface dont une partie au moins 
represente la surface d'un obstacle recevant une onde 
principale et emettant, en reponse, une onde secondaire, 
et l'on attribue a chaque echantillon de surface au moins 
une source emettant une onde elementaire representant une 
contribution a ladite onde secondaire, 

b) on forme un systfeme matriciel comport ant : 

- une matrice d' interaction, inversible, appliguee a 
une region donnee de 1 ' espace et comport ant un nombre 
de colonnes correspondant a un nombre total de 
sources, 

- une premiere matrice colonne dont chaque coefficient 
est associe a une source et caracterise l'onde 
elementaire qu'elle emet, 

- et une seconde matrice colonne, obtenue par une 
multiplication de la premiere matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d'une grandeur physique representative de 
l'onde emise par 1' ensemble des sources en ladite 
region donnee, 

O pour estimer les coefficients de la premiere matrice 
colonne, on affecte des valeurs de grandeur physique 
choisies a des points predetermines, associes chacun a un 
echantillon de surface, lesdites valeurs choisies etant 
rangees dans la seconde matrice colonne, et l'on multiplie 
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cette seconde matrice colonne par 1 ■ inverse de la matrice 
d« interaction appliquee auxdits points predetermines, 
d) pour evalue:r la^'te grandeur physique representant 
l'onde emise par l 1 ensemble des sources dans line region 
5 donnee de l'espace tridimensional, on applique la matrice 
d' interaction a ladite region donnee et on multiplie cette 
matrice d' interaction par la premiere matrice colonne 
comportant les coefficients estimes & l'etape c) . 

10 Ainsi, selon 1'un des avantages que procure la presente 
invention, l/etape de maillage a) ne concerne qu'une ou 
plusieurs surfaces, tandis que le proc6de de modelisation 
de type "par elements finis" ngcessite un pavage de tout 
l'espace avoisinant 1' obstacle, ce qui permet de reduire, 

15 dans la mise en oeuvre du proced§ selon 1" invention, les 
ressources mgmoires et les temps de calcul n^cessaires. 

Le precede selon l 1 invention s' applique aussi bien a une 
oride principale Smise par une source lointaine qu'& une 
20 onde principale emise en champ proche. 

Avantageusement, pour ^valuer une grandeur . physique 
representative d'une interaction entre un element 
rayonnant une onde principale en champ proche et un 

25 obstacle recevant cette onde principale, 

- a l'etape a), on determine en outre, par maillage, une 
plurality d'echantillons de surface representant ensemble 
une surface active de 1« element rayonnant l'onde 
principale, et I'on attribue & chaque echantillon de la 

30 surface active au moins une source Smettant une onde 
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§l€mentaire repr§sentant une contribution a ladite onde 
principale, 

- on applique les ecapes b) , c) et d) aux Schantillons de 
la surface active, et 
5 on evalue ladite grandeur physique repr^sentant 

1 1 interaction entre l'61£ment rayonnant et 1 1 obstacle dans 
une region donn£e de l'espace tridimensionnel, en prenant 
en compte la contribution, en ladite region donn#e, de 
I'onde principale 6rnise par I 1 ensemble des sources de la 
10 surface active et la contribution de l'onde secondaire 
gmise par 1 1 ensemble des sources de la surface de 
1 'obstacle. 

On entend par les termes "element rayonnant" aussi bien un 
15 emetteur de l'onde principale, tel qu f un g£nerateur 
d'ondes, qu'un recepteur de l'onde principale, tel qu'un 
capteur de cette onde. 

Dans un premier mode de realisation, la grandeur physique 
20 & evaluer est une grandeur scalaire et, & 1 1 £tape a), on 
attribue une source unique a chaque gchantillon de 
surface. 

Dans un second mode de realisation, la grandeur physique & 
25 §valuer est une grandeur vectorielle exprimSe par ses 
trois coordonnSes dans 1 1 espace tridimensional, et on 
attribue, a 1 1 etape a), trois sources a chaque echantillon 
de surface. 

30 Dans une realisation avantageuse, pour estimer, a 1' etape 
d) , la contribution de l ! onde secondaire dans la region 
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donnee de l'espace, les valeurs de grandeur physique 
choisies a l'etape c) sont fohction d'un coefficient 
predetermine de reflexion et/ou de transmission de l'onde 
principale par chaque echantillon de surface de 
1 'obstacle. 

Ainsi, on comprendra que l'onde secondaire peut aussi bien 
correspondre a une reflexion de l'onde principale, qu'a 
une transmission de l'onde principale, ou encore a une 
diffraction de l'onde principale. Dans cette realisation 
avantageuse, l'etape c) correspond finalement a une 
determination des conditions aux limites a la surface de 
1 'obstacle, en tant qu ' interface entre deux milieux 
distincts notamment dans une hSterostructure . 



Par ailleurs, pour un contrdle non destructif d'une cible 
formant un obstacle d'une onde principale, on attribue un 
coefficient de reflexion ou de transmission choisi a tous 
les points predetermines de la surface de la cible, et on 
compare une simulation obtenue par la raise en ceuvre du 
procSde au sens de 1' invention avec une mesure 
experimentale. Ainsi, les points de la surface de la cible 
qui, dans la mesure experimentale, ne verifient pas la 
simulation correspondent a des inhomogeneites ou a des 
25 impuretes sur la surface de la cible. 

Dans une autre approche, on connait les proprietes 
globales de 1' obstacle, notamment en transmission et/ou en 
reflexion. Par la mise en oeuvre du proc£d£ de 1' invention, 
30 on optimise alors la position dans l'espace d'un capteur 
ou m§me la forme de ce capteur, pour une application a un 
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controle non destructif , ce capteur «taat destine a 
analyser une cible formant un obstacle de l'onde 
principale. 

A cet effet, dans une realisation avantageuse, on compare 
une plurality de valeurs de la grandeur physique estimee a 
l'etape d) du precede au sens de 1- invention, obtenues 
pour une plurality de regions de 1'espace, pour 
selectionner une region candidate pour la disposition d'un 
element rayonnant destine a interagir avec 1- obstacle. 

Comme indigue ci-avant, on entend par les termes "element 
rayonnant" aussi bien un capteur qu'un generateur de 
l'onde. On comprendra ainsi que 1 - optimisation de la 
position de I- element rayonnant pent s'appliquer aussi a 
1' optimisation de la disposition ou de la forme d'un 
generateur d'ondes. Par exemple, l a presente invention 
trouve encore une application avantageuse a la disposition 
de haut-parleurs dans un volume ferme, delimite par des 
obstacles, comme par exemple l'habitacle d'un vehicule 
automobile. 

D'autres caracteristiques et avantages de 1' invention 
apparaitront a l'examen de la description detaillee ci- 
apres, et des dessins annexes sur lesquels : 
- la figure 1A represente schematiquement les surfaces 
respectives d'un element rayonnant ER emettant une onde 
et d'un obstacle DBS recevant cette onde, maillees en 
vue d'evaluer une grandeur scalaire representative de 
l'onde en un point M de l'espace tridimensionnel ,- 
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- la figure IB represente en detail un echantillon de 
surface dS± correspondant a une' maille de la figure 1A, 
ainsi qu'une source Si associee a 1 'echantillon de 
surface dSi ; 

5 - la figure 2A represente schematiquement les surfaces 
respect ives d'un element rayonnant ER emettant une onde 
et d'un obstacle OBS recevant cette onde, maillees en 
vue d'evaluer une grandeur vectorielle representative 
de l'onde en un point M de l'espace tridimensionnel ; 
10 - la figure 2B represente en detail un echantillon de 
surface dS £ correspondant a une maille de la figure 2A, 
ainsi que trois sources associees SAi, SBi et SCi ; 

- la figure 2C represente, en vue.de face, une surface 
maillee dont chaque 6chantillon de surface comporte 

15 trois sources SAi, SBi et SCi, pour 1 -estimation d'une 

grandeur vectorielle ; 

- la figure 3A represente, a titre illustratif, les 
armatures d'un condensateur , de potent iels electriques 
respectifs VI et V2, pour 1' estimation d'un potentiel 

20 electrique au point M de l'espace tridimensionnel, a 

chaque echantillon de surface dS d de la figure 3A etant 
associee une unique source Si ; 

- la figure 3B represente, a titre illustratif, les 
armatures d'un condensateur, de potentiels electriques 

25 respectifs VI et V2, pour 1 "estimation d'un champ 

electrique E(M) , au point M de l'espace 
tridimensionnel, a chaque ichantillon de surface dSi de 
la figure 3B etant associee trois sources SAi, SBi et 
SCi ; 
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la fxgure 4A represente, comme les figures 1A et 2A, 
une interaction entre un element rayonnant ER et un 
.obstacle OBS, pour ^ luer ^ grandeur phys±que 
(scalaire ou vectorielle) en un point M dans une 
portion de 1 ■ espace delimits par la surface de 
1' element rayonnant et la surface de 1- obstacle, ce 
point de I- espace M recevant a la fois 1'onde emise par 
1- element rayonnant et 1-onde reflechie par 
1 ' obstacle ; 

10 - la figure 4B, complement aire de la figure 4A 
represente une transmission par 1 'obstacle OBS de 
1'onde emise par 1- element rayonnant ER, en un point M 
d-un demi-espace delimits par le plan que forme la 
surface de 1- obstacle OBS ; 

15 - la figure 5A represente schematiquement un obstacle 
OBS, de dimensions finies, avec des sources associees 
aux echantillons de surface agencees pour estimer une 
grandeur representative d'une reflexion de 1-onde sur 
1' obstacle ; 

20 - la figure 5B, en complement de la figure 5A, represente 
schematiquement un obstacle OBS, de dimensions finies 
avec les sources associees aux echantillons de surface 
dxsposees pour estimer une grandeur representative de 
la transmission de l'onde par 1 'obstacle ; 

25 - la figure 5C represente une simulation d'une onde 
ultrasonore emise par un element rayonnant ER et se 
propageant vers un obstacle OBS ; 

la figure 6 represente un maillage de milieux 
tridimensionals, pour 1 -application d-un precede de 
calcul par "events «„,.., au sens de ^ 
technique ; 
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la figure 7A represente en detail un element de surface 
et un point d' observation M dont les positions 
relatives sont reperees par un angle 6 ; et 
la figure 7B represente schematiquement une surface a 
5 mailler de forme complexe, en particulier avec un point 

d- observation M se situant dans une zone d- ombre par 
rapport a certaines sources de la surface. 

On se refere tout d'abord a la figure 1A, sur laquelle la 
10 surface d'un obstacle OBS, recevant une onde, est maillee 
selon une plurality d' echantillons de surface d Sl a dS 4 , 
conformement a l'etape a) precitee. 

En se referant a la figure IB, a chaque echantillon de 
15 surface d Si est associe un hemisphere HEMi, tangent a 
1' Echantillon de surface dS* en un point de contact P A . 
Preferentiellement, ce point de contact P, correspond au 
sommet de 1- hemisphere HEMi . Pour !• estimation d-une 
grandeur physique scalaire au point M (tels qu'un 
20 potentiel electros tat ique, une pression acoustigue ou 
autre), une source unique Sl est associ6e a 1 ■ echantillon 
de surface d Si . Comme on le verra plus loin, dans le cas 
de 1 'estimation d'une grandeur vectorielle dans un point 
de 1'espace M, on affectera P lut6t trois sources a chaque 
25 Echantillon de surface dSi. 

Preferentiellement, 1 -hemisphere HE Mi est construit comme 
decrit ci-apres. Pendant l'etape de maillage a) precitee, 
on evalue, d'une part, la surface de 1' obstacle OBS et, 
30 d' autre part, on choisit un nombre d ' echantillons de 
surface dSi selon la precision souhaitee de 1' estimation 
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de la grandeur physique au point M. Ainsi, la surface d'un 
echantillon dSi est donn€e par S Q /N, oil S 0 correspond a la 
surface totale de 1 1 obstacle et N correspond au nombre 
choisi d'echantillons de surface dSi. 

5 

L ' hemisphere HEMi est de meme surface que 1 1 echantillon 
dSi. Ainsi, le rayon R t de 1 ■ hemisphere se deduit de 
.1 'expression : 

1 N 

Chaque maille que represente un echantillon de surface dS± 
presente, dans l'exemple decrit, une forme de 
paralieiogramme, de centre P ± correspondant au point 
d 1 intersection des diagonales de ce paralieiogramme . 
15 L 1 hemisphere HEMi est tangent k 1 ' echantillon de surface 
dSi en ce point Pi. Bien entendu, les mailles peuvent §tre 
de forme differente, triangulaire ou autre. On indique de 
fagon gen^rale que le point Pi correspond au barycentre de 
la maille . 

20 

On definit ainsi la position de la source Si (situ^e au 
centre de 1' hemisphere HEMi) • La distance separant la 
source Si du point de contact Pi correspond au rayon R A de 
1 'hemisphere HEMi et la droite qui passe par les points Pi 
25 et Si est orthogonale a la maille dSi. 

Dans I'exemple represente sur la figure 1A, on maille en 
outre la surface d'un element rayonnant ER, correspondant 
par exemple a un generateur d'ondes. A chaque maille de la 
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surface de 1> element rayonnant est associ<§ un echantillon 
de surface dS' if comme on le verra plus loin. 

Le systeme matriciel que l'on met en forme a 1 ' etape b) 
precitee correspond a : 



^(M!) > 






V(M 2 ) 


= F x 


V2 


.V(M N ), 







[1] 



ou 



10 



15 



20 
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les coefficients v., (avec j=l,2,..,N) de la premiere 
matrice colonne correspondent a des valeurs associees 
chacune a une source, telle qu'une charge electrique 
(dans le cas de 1 -estimation d'un potentiel 
Electrique) , ou encore a un flux magnetique (dans le 
cas de 1' estimation d'un potentiel magnetique), ou 
encore a une vitesse sonore (dans le cas de 
1' estimation d>une pression acoustique liee a la 
propagation d'une onde sonore) ; 

les coefficients V( Mi ) (avec i=l,2,.„,N) de la seconde 
matrice colonne correspondent chacun a une valeur de la 
grandeur physique (un potentiel Slectrique ou 
magnetique ou une pression) a estimer en un point M± de 
1 1 espace ; 

la matrice d- interaction P comporte des coefficients 
Ci,j dont 1' expression generale est donnee par : 



Gj=f{M,Sj) 



[2] 
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On comprendra ainsi que les coefficients de la matrice F 
sont des coefficients d' interaction qui dependent de la 
distance separant chaque point de l'espace M ± d'une source 
Sj associee a une maille dSj . 

Dans le cas de la propagation d'une onde Slectrique, les 
coefficients ci.-j, Vj et V(M ± ) , respectivement de la 
matrice d' interaction de la premiere et de la seconde 
matrice colonne, sont donnas par : 



r 1 
c i,j = 



2xe 0 MjSj 
V(Mi) = Ui 

oil : 

- s 0 est une constante diglectrique, 

- MjSj est une distance mesuree en valeur algehrique, 

- qj correspond a une charge electrique caracterisant une 
source Sj, et 

- Ui correspond a un potentiel Electrique au point Mi. 

Dans le cas de la propagation d'une onde magnetique, 
1' expression de ces coefficients est la suivante : 



C 

1,3 2n Mo MjSj 



V j =< Pj [4] 



V(Mi) = ^ 
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oil : 

- fi a correspond a la pexmeabiiice magn£tique du milieu ou 
se situe le point Mi, 

5 - <pj correspond au flux magnetique associ€ a la source 
Sj ; 

- &x correspond au potentiel magngtigue au point Mi. 

Dans le cadre de la propagation d'une onde ultrasonore, 
10 ces coefficients sont donnes par : 



_ leap expl 



(ikJvljSi) 

. ua ; 

2n M^j J 



Vj=Vj 



t5] 



V(Mi)= Pi 

dans laquelle : 
15 - i 2 =-l. 

- co est la pulsation de l'onde sonore ; 

- p est la density du milieu dans lequel se situe le 
point Mi ; 

- le vecteur Vj correspond a la vitesse sonore issue de 
20 la source Sj ; 

k correspond au vecteur d'onde de l f onde sonore ; et 

- Pi correspond & la pression acoustique g€n#ree par la 
propagation de l'onde ultrasonore au point Mi. 
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Dans 1 'expression des coefficients e lfJf le terme dSj 
correspond a la surface de 1 - echantillon associe a la 
.source S 3 . De preference, le maillage d'une surface au 
sens de l'etape a) du precede selon 1- invention, est 
choisi de sorte que chaque maille comporte une m§me 
surface dS=dS 1 =dS 2 =...=dS j . 

On remarque en particulier dans 1- expression des 
coefficients c lfi qu'ils dependent du produit scalaire 
entre le vecteur d'onde et le vecteur . Ainsi, pour 

des ondes ultrasonores, on tient compte d'un dephasage 
entre les chemins qui lient chaque source S d a un point de 
l'espace tridimensional M, ce dephasage etant du a une 
difference de marche entre les rayons partant de chaque 
source et arrivant au point M (comme le tnontre la figure 
4A). En particulier, on comprendra que 1- angle d- incidence 
d'un tel rayon est pris en compte dans 1 • expression des 
coefficients de la matrice d 'interaction P. 

Bien entendu, dans le cadre de la propagation d-une onde 
elect romagnetique de fr6quence elevee, done de longueur 
d-onde courte, qui differe du cadre electrostatique ou 
magnetostatigue ci-avant, on tient compte du terme de 
propagation exp^.r) dans 1- expression de la matrice 
d- interaction, par rapport a la geometrie du probleme a 
resoudre, comme dans le cadre de la propagation d'une onde 
ultrasonore ci-avant (relation [5] ) . 

Ainsi, le systeme matriciel de 1- equation [1] permet 
d'estimer, a partir d'une matrice d- interaction F et d'un 
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vecteur comportant les valeurs 
Sj, les coefficients d'un 



associ6es aux sources 
vecteur (matrice colonne) 
comportant les valeurs de grandeur physique V(M A ) aux 
points de l'espace M±. 



Pour determiner les valeurs des sources v jf on applique le 
systdme matriciel suivant : 



10 



15 



20 



les 



r V(Pi) > 








V(P 2 ) 


- F x 




[6] 










coefficients 


de 


la matrice d ' interaction P 



s'expriment par avec Gj=f[BSj) 
les indices i et j correspondent respect ivement a la 
i ftme ligne et la j*" e colonne de la matrice d' interaction 
F. Cette matrice d« interaction comporte, pour la 
determination des valeurs associees aux sources Vj , N 
lignes et N colonnes, en rappelant que N est le nombre 
total de mailles a la surface de 1' obstacle ; 
les points Pi correspondent au sommet des hemispheres 
HEMi de la figure IB. 



25 



La mise en cfiuvre de l'etape c) du precede au sens de la 
presente invention correspond a calculer une condition aux 
limites pour les points Pi , de proprietes connues, comme 
on le verra plus loin. 
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Le systSme matriciel de liquation [6] devient alors : 











= F _1 x 


V(P 2 ) 






,V(P N ), 



[7] 



5 oil : 

- F" 1 correspond a 1 ■ inverse de la matrice d' interaction 
F ; et 

- les valeurs V(Pi) sont predetermines, en fonction des 
conditions aux limites precitges. 

10 

On determine ainsi les valeurs de source v j . 

A partir de l 1 estimation de ces valeurs de source v-j, on 

peut calculer la grandeur physique scalaire en un point M 

15 quelconque de 1 ' espace tridimensional, a partir de la 
relation : 

V(M)=^MSj)vj [8 ] 

20 Pour obtenir cette expression de la grandeur scalaire 
V(M), la matrice d' interaction F peut ne comporter qu'une 
ligne de coefficients Cj, avec : 

Cj = £(MSj) , 
mais comporte toujours N colonnes. 

25 

En se r€f grant k nouveau §l la figure 1A, on comprendra que 
la surface de 1" obstacle OBS recevant I'onde qu'€met 
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1 Element rayonnant ER agit, elle-meme, comme une surface 
active re-emettant une onde secondaire (par exemple par 
reflexion) . chaque source S d represente une contribution a 
1' emission de cette onde secondaire. 



10 



En outre, pour tenir compte a la fois de la presence de 
l'onde principale et de la presence de 1'onde secondaire 
au point M, on estime la contribution de 1-onde principale 
et la contribution de l'onde secondaire au point M par le 
systeme matriciel : 



V(M) = Fx 



fv, > 




fv'i ' 










+ Fx 











[9] 



15 



20 



OU 



F' est la matrice d • interaction entre la surface de 
1' element rayonnant et le point ? M ; 

V ' 3 (avec 3=1,2,3 N-) est la valeur des sources 

attributes a chaque e chant il Ion de surface dS'j de 
1' element rayonnant, N- etant le nombre total de 
mailles choisi pour la surface active de 1- element 
rayonnant ER. 



25 



Les coefficients de la matrice P- sont encore fonction de 
la distance MS> it ou S'j sont les sources affectees a 
chaque echantillon de surface dS-j de !• element rayonnant. 



Selon une caracteristique avantageuse, les valeurs des 
sources de 1 • obstacle Vj sont determines en fonction des 
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valeurs des sources de 1' element rayonnant v' j# lesquelles 
sont elles-memes calculees comme on le verra plus loin en 
reference aux figures 4A, 4B, 5A et 53. 

5 On se refere maintenant a la figure 2A, dans laquelle on 
affecte trois sources a chaque echantillon de surface dS i# 
en vue d'estimer une grandeur physique vectorielle V(M) , 
en un point M de l'espace tridimensional . 

10 En effet, on comprendra que pour estimer la grandeur 
vectorielle, par ses trois coordonnees dans l'espace x, y 
et z, le nombre d' equations a resoudre par rapport au 
systeme matriciel precedent doit etre multiplie par trois. 
Ainsi, la matrice F" 1 de la relation [7] doit comporter 
trois fois plus de lignes que precedemment. La matrice 
d' interaction P doit, elle-meme, comporter trois fois plus 
de colonnes que precedemment et, a cet effet, on prevoit 
avantageusement trois sources par maille lorsqu'il s'agit 
de determiner les coordonnees dans l'espace 
tridimensionnel d ' un vecteur V(M) . 



15 



20 



25 
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En se referant a la figure 2B, les trois sources SAt, SBi, 
SCi, attributes a un echantillon de surface dSi sont de 
positions respect ives determinees comme on l'indique ci- 
apres. Telles que representees sur la figure 2B, les trois 
sources SA A , SBi, Sd sont coplanaires et le plan 
comport ant ces trois sources comport e en outre la base de 
1 'hemisphere HEMi . L- hemisphere HEM* est construit comme 
indiqug ci-avant (de merae surface que la surface de la 
maille), avec toutefois le centre de 1' hemisphere qui 
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correspond ici au barycentre des trois sources SAi, SBi et 
SC ± . 

On entend par "centre de 1 'hemisphere" le centre du disque 
qui constitue la base de 1 'hemisphere . 

Les trois sources qui sont attributes a 1 ' echantillon de 
surface dSi sont disposees aux sommets d'un triangle 
equilateral dont le barycentre Gi correspond au centre de 
1 'hemisphere. De preference, chaque source SAi, SBi et SC ± 
est disposee sur le milieu d'un rayon R ± del' hemisphere. 
Ainsi, les droites qui relient le barycentre Gi a chaque 
source sont ecartees angulairement de 120°. 

En se referant a la figure 2C, 1 ' orientation angulaire des 
triangles que forment les triplets de sources est choisie 
aleatoirement, d'un echantillon de surface a 1' autre. 
Avantageusement, on evite ainsi des artefacts de sur- 
periodicite, dans 1' estimation de la grandeur vectorielle 
au point M, qui pourraient resulter du choix d'une meme 
orientation angulaire de ces triangles. 



En reference avec les differents types d'ondes indiques 

precedemment , la grandeur vectorielle V(M) a estimer peut 
25 etre : 

- un champ electrique, dans le cadre de la propagation 
d'une onde Electrique ; 

- un champ magnet ique, dans le cadre de la propagation 
d'une onde magnetique ; et 

- une Vitesse du son au point M, dans le cadre de la 
propagation d'ondes ultrasonores . 
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Pour determiner les valeurs associ€es a chaque source SA±, 
SB it SCi, on met en forme le systeme matriciel selon la 
relation suivante : 



r V x (P!)> 
V X (P 2 ) 




"vA x ^ 
VA2 


V X (P N ) 
V y (Pi) 

Vy(P 2 ) 


= Fy X 


vA N 

vBj 

vB 2 


Vy(P N ) 

V z (Pi) 
V 2 (P 2 ) 




vB N . 

vCj 

vC 2 


^(Pn); 







[11] 



10 



On remargue, en particulier, que la matrice d' interaction 
Fy est de dimensions 3Nx3N, 01a N est le nombre total 
d'6chantillons de surface. La matrice d » interaction 
s'exprime ici par la relation : 



F V 



N 



N 



N 



NRftj) C|(i,j) C£(i,j) 
N^A^J) Cj0.j) Cl(i,j) 
N{= A (i.j) Cl(i,j) C£(i,j) 



[12] 
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Les coefficients de cette matrice s 1 expriment par : 

cZ(ij)=f u [d(p i9 s*jj\ tl31 

5 avec a = A, B, C 

i - 1, 2, N 
j - 1, 2, N 
u = x, y, z. 

10 et sont encore fonction d*une distance sgparant le point 
de contact Pi de l'une des sources Soj (a = A, B ou C) 
d»un triplet associe a un echantillon de surface dSj . 

Par inversion de la matrice d 1 interaction , Fy , on 

15 determine ainsi les valeurs vcrj associees & chaque source 

Saj, en appliquant des conditions aux limites sur les 

valeurs du vecteur V aux points Pi. Ces conditions aux 

limites imposent une valeur du vecteur V , selon ses trois 
coordonn§es V x (Pi) , V y (Pi) et V z (Pi). 

20 

Une fois que ces valeurs de sources vcrj sont ainsi 
d§termin€es, on calcule facilement I 1 expression du vecteur 
V en un point quelconque M de l'espace, par la relation : 

25 V(M) = V X (M) x + V y (M) y + V 2 (M) z [14] 

cr=A,B,C 
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J=l,..,N 
<y=A,B,C 



10 



7=1 N 

a=A,B,C 

x, y et z correspondent a des vecteurs unitaires portes 
par les axes x, y et z de 1 > espace tridimensionnel . 

Ainsi, la matrice d • interaction Fy , lorsqu'elle est 
appliquee a un point guelconque M de 1 -espace, ne comporte 
finalement que trois lignes associees chacune a une 
coordonnee de 1 ' espace x, y ou z. 

Pour differents types d'ondes, les valeurs des sources vcr-j 
sont, comme prScedemment , une charge Slectrique pour une 
onde electrique, un flux magnet igue pour une onde 
magnetique, une vitesse du son pour une onde ultrasonore. 

Plus precisement, les coefficients de la matrice 
d' interaction F y se determined a partir des relations 
[3] , [4] et [5] precedentes en precisant toutefois que : 

V(M) = -i^d[V(M)] [15] 

V(M) etant la grandeur scalaire calculee precedemment par 
25 1 'equation [8] . 
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Ainsi, pour 1 1 estimation d'une grandeur vectorielle V au 
point M et pour les types d'onde cites ci-avant a titre 
d'exemple (electrique, magnetique et ultrasonore) , les 
coefficients de la matrice d' interaction Fy sont 
5 inversement proportionnels au carre d f une distance * 
sEparant chaque source du point M, alors que pour 
1 1 estimation d'une grandeur scalaire V en un point M de 
l'espace, les coefficients de la matrice d 1 interaction F 
sont simplement inversement proportionnels a cette 

10 distance. Chaque distance implique l'une des sources d'un 
triplet d'un echantillon de surface et un point M de 
l'espace. La matrice d' interaction comporte alors 3N 

colonnes lorsqu'il s'agit de prendre trois sources par 
echantillon de surface, alors que la matrice d' interaction 

15 F pour I 1 estimation de la grandeur scalaire ne comportait 
que N colonnes puisque seule une source par Echantillon de 
surface etait nEcessaire. • 

Plus generalement, on attribue une source par echantillon 
20 lorsque l'on connait des conditions aux limites pour une 
grandeur scalaire et trois sources par Echantillon lorsque 
l'on connait plut6t des conditions aux limites pour une 
grandeur vectorielle. 

25 On se r£fere maintenant a la figure 3A pour decrire, a 
titre illustratif, une application du procEdE selon 
1' invent ion a l f estimation d'un potentiel electrique en un 
point M de I'espace tridimensional, situe entre deux 
armatures d'un condensateur. Les armatures de ce 

30 condensateur sont portEes a des potentiels respectifs VI 
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et V2. La mise en oeuvre de l'etape a) consiste d'abord a 
mailler les surfaces respectives des deux armatures, Dans 
l'exemple represent! aur la figuxe 3A, on n'a repffesente 
que deux mailles pour chaque armature, simplement a titre 
illustratif. 



10 



L» application de l'etape b) consiste a mettre en forme le 
systeme matriciel impliquant la matrice d 1 interaction F et 
le vecteur colonne comportant les valeurs des sources Si & 
S 4 . La multiplication de ces deux matrices permet 
d f obtenir un vecteur colonne comportant les valeurs du 
potentiel en un ou plusieurs points M de l'espace. 



15 



La mise en oeuvre de l f §tape c) du procede selon 
l 1 invention consiste £ appliquer le systeme matriciel aux 
points de contact des hemispheres Pi a P 4 , de chaque 
echantillon de surface dS x & dS 4 . 11 en r£sulte la 
relation suivante : 



20 



v(p 2 ) 
v(p 3 ) 
Wp 4 )J 



2k£q 



PA 

1 


P1S2 

1 


P1S3 

1 


P1S4 

1 ! 


1 


P 2 S 2 

1 


P 2 S 3 

1 


P 2 S 4 

1 


P3S1 

1 


P3S21 
1 


P3S3 
1 


P3S4 

1 


P4S1 


P4S2 


P4S3 


P4S4; 



V3 



[16] 



avec 



V(Pi) = V(P 2 ) = Vi 
V(P 3 ) = V(P 4 ) = V 2 
Vi = q x , v 2 = q 2 , v 3 



q3/ v 4 



q4 
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Ici, la condition aux limites impose que la valeur du 
potentiel aux points de contact P 1 et P 2 correspond au 
potentiel VI de la premiere armature. De meme, le 
potentiel Slectrique aux points de contact P 3 et P 4 
correspond au potentiel Slectrique de la seconde armature 
V2. Par inversion de la matrice d- interaction appliquee au 
point de contact P if on determine les valeurs des sources 
v ± qui correspondent, comme exprime dans la relation [16] , 
a des charges electriques q A . 

Les coefficients de la matrice d ' interaction l - SO nt 

2^ 0 PiSj 

parfaitement connus, puisque les positions des sources Sj 
et les positions des points de contact P 4 sont determines 
au prealable, comme represents sur la figure IB. 

L ' expression du potentiel electrique V(M) au point M entre 
les deux armatures est finalement donnSe par 
1' expression : 

V(M) = -J-f-Sl_ + -52_ + _a3_ + _a4_ > | ri7 , 
2ne 0 { MS, MS 2 MS 3 MS 4 J U7j 

On se rSfere maintenant a la figure 3B sur laquelle on a 
represents les memes armatures que sur la figure 3A, avec 
sensiblement un m§me maillage, mais dans le but, ici, 
d'estimer une grandeur vectorielle correspondant au champ 
electrique E(M) , au point M de l'espace tridimensionnel . 
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On peut appliguer les relations [11] a [15] pour estimer 
la valeur du champ electrique au point M, avec, dans la 
relation [13] : 



f 



2ns Q d 2 {P u S<jj)^ 



[18] 



avec 

c = A, B, c 

u = x, y, z 

i - 1, 2, 3, 4 

10 j = 1, 2, 3, 4. 



Toutefois, les valeurs du champ Electrique au point de 
contact Pi restent a determiner dans la relation [11] . 

On introduit alors une loi generale prEdeterminee du 
comport ement du champ (en reflexion, en transmission, ou 
autre) au niveau de la surface de 1' obstacle (des 
armatures dans l'exemple du condensateur precite) , pour 
connaitre les valeurs des sources Vaj . 



Par exemple, si l'onde electrique est totalement reflechie 
par la surface d'un obstacle (par exemple l'une des deux 
armatures) , le champ Electrique en un point de contact Pi 
est normal & la surface dSi et ses composantes E x et Ey 
25 sont nulles. A titre illustratif, si la surface de 
1' armature n'etait representee que par un seul echantillon 
de surface a trois sources, les valeurs de ses sources vA, 
vB et vC seraient toutes egales ent're elles a une m§me 
valeur +q. 
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Au contraire, si le coefficient de reflexion est 
pratiquement nul a la surface dS x , la composante du champ 
Slectrique E z au point Pi est nulle, ce qui correspond 
bien au cas oil le champ est sensiblement tangent a la 
surface dSi. Ainsi, a titre illustratif, si la surface de 
1' armature n'etait representee que par un seul echantillon 
de surface a trois sources, les valeurs de ses sources vA, 
vB et vc seraient respect ivement, par exemple, +q, +q et 
-2q. Par exemple, dans le cadre de la propagation d'une 
onde magnet ique, si la surface d'un capteur a courants de 
Foucault (avec une composante normale du champ magnetique 
nulle) etait representee par un unique echantillon de 
surface, les flux magnetiques des trois sources associees 
a cet echantillon de surface seraient +9, +<p et -2<p. 

On comprend ainsi qu'avec les trois sources par 
echantillon dSi, on peut definir, par exemple en fonction 
de la ponderation de chaque source, une orientation 
quelconque du champ a la surface de 1' obstacle. 



Bien entendu, cette demarche suppose que le coefficient de 
reflexion R d'un obstacle soit connu auparavant. En 
particulier, il peut etre avantageux de comparer une 
simulation et une mesure experimental pour detecter, a la 
surface d'un obstacle, des inhomogeneites ou des impuretes 
qui correspondent a des points de la surface de cet 
obstacle qui ne verifient pas les valeurs du coefficient 
de reflexion R imposees a chaque point predetermine Pi de 
30 1' obstacle. 
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On peut af fecter ainsi une valeur pred€termin£e . du 
coefficient de reflexion & chaque point Pi de la surface 
de 1' obstacle. A cet effet, on incroduit une matrice R qui 
est representative du coefficient de reflexion en chaque 

5 point Pi, Pour une interaction entre un Pigment rayonnant 
et un obstacle, on peut ainsi exprimer le systeme 
matriciel de la relation [9] autrement, c'est-a-dire en 
dormant une expression -unique de toutes les sources du 
systeme (k la fois de l 1 obstacle et de 1 "filament 

10 rayonnant), comme indiqu6 ci-apr£s. 

Dans ce qui suit, on indique que : 

- F(P) est la matrice d ! interaction de 1' obstacle OBS 
appliquee aux points Pi de la surface de l 1 obstacle 

15 OBS ; 

- F(P') est la matrice d ! interaction de 1" obstacle OBS 
appliquee aux points P'i de la surface de I 1 element 
rayonnant ER; 

- F 1 (P) correspond & la matrice d 1 interaction de 
20 1" element rayonnant ER appliquee aux points Pi de la 

surface de 1 1 obstacle OBS ; 
• F' (P f ) correspond k la matrice d 1 interaction de 
1 1 Pigment rayonnant ER appliquee aux points P± de la 
surface de l'§16ment rayonnant ER ; 
25 - v* correspond au vecteur colonne comportant les valeurs 
des sources S 1 i de 1 1 £l€ment rayonnant ER ; et 
v correspond au vecteur colonne comportant les valeurs 
des sources Si de l 1 obstacle OBS . 
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Sur 1* obstacle, la contribution de l'onde emise par 
1 Element rayonnant ER s ' exprime par : 

V'(P) = F(P).v' [19] 

La contribution de l'onde secondaire renvoyee par 
1' obstacle OBS s' exprime, par definition, par la 
relation : 



10 V (P) = F(P).v [2 0] 

Or, dans 1 1 exemple reprgsente sur la figure 4A, l'onde 
secondaire correspond simplement £ une reflexion de l'onde 
principale. Ce qui s 1 exprime par la relation : 



15 



20 



V(P) = RV'(P) [21] 

oil R correspond & une matrice de reflexion dont chaque 
coefficient reprisente la contribution k 1' Emission, par 
reflexion, de l'onde secondaire, par chaque source Si (ou 
So it dans le cadre d'une estimation d'une grandeur 
vectorielle) de 1' obstacle OBS . 



Des trois relations [19] , [20] et [21] , on d^duit 
25 l 1 expression du vecteur colonne v comportant les valeurs 
des sources de 1 'obstacle, a partir du vecteur colonne v 
comportant les valeurs des sources de 1 1 Pigment rayonnant, 
par la relation : 



30 



v = [F(P)]' 1 .R .[P i (P)].v i 



[22] 
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Par ailleurs, pour une estimation fine des grandeurs 
scalaires ou vectorielles au point M, nocamment pour tenir 
compte de reflexions multiples, il est avantageux de tenir 

5 compte de la contribution du rayonnement par 1 ! obstacle, 
au niveau de la surface de l'^l^ment rayonnant ER. A cet 
effet, on tient compte, dans 1 ' estimation des conditions 
aux limites a la surface de 1' Element rayonnant ER (aux 
points P ! i) de la contribution du rayonnement des sources 

10 S § i de l 1 element rayonnant et de la contribution de 
1' Emission de l'onde secondaire par les sources Si de 
I 1 obstacle, par la relation : 

VtCPO^FCPOv + FCF)^ [23] 

15 

On peut ainsi ajuster, grSce. k la relation [23], les 
valeurs de sources S'i de l 1 Element rayonnant ER, en 
tenant compte de la reflexion de 1" obstacle OBS, selon la 
relation suivante : 

20 

V T (P') = \f(P') .[F(P) ]~\ .R .[F'(P) ] + F'(P') }.v' [10] 

Ainsi, on impose simplement des conditions aux limites 
pour 1' Element rayonnant, pour en dfiduire les valeurs des 
25 sources v'i. En pratique, on procedera pref €rentiellement 
comme suit : 

- apres maillage des surfaces, on determine la position 
des points P A et P ! t et des sources S± et S'i; 

- en fonction du type d f onde en jeu, on determine les 
30 coefficients des matrices F(P), F 1 (P) , F(P f ) et F'tP 1 ); 
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- en fonction d'une loi de reflexion de 1 -obstacle, on 
determine les coefficients de la matrice de reflexion 
comme dans l'exemple donne plus loin pour une onde 
ultrasonore; 

- en fonction de conditions aux limites sur 1' element 
rayonnant (dont le comportement est generalement connu 
pour un probleme donne), on determine les valeurs du 
vecteur V«F) aux points P'i de la surface de 1- element 
rayonnant et on en deduit les valeurs des sources S'i de 
1- element rayonnant par inversion de la relation [10]; 

- on en deduit aussi les valeurs des sources S 4 de 
1' obstacle par application de la relation [22]; 

- une fois les valeurs de toutes les sources S ± et S'i 
determines, on peut appliguer le systeme matriciel donne 

15 par la relation [9] a tout point M de l'espace, en 
appliquant a ce point M les matrices d- interaction F et F' 
(impliquant la position du point M et les positions des 
sources respectives Si et S'i). 



10 



20 



25 



En se referant a nouveau a la figure 4A, on considere que 
1' obstacle OBS represente simplement une interface entre 
deux milieux Ml et M2, formant ainsi un dioptre qui peut 
etre plan, tel que represente dans l'exemple de la figure 
4A, mais aussi incurve ou de forme generale quelconque. 
Les coefficients de reflexion Rl associes a chaque point 
Pi dependent, dans le cadre de la propagation d'une onde 
ultrasonore ou electromagnetique de haute frequence, de 
1' angle d- incidence 0 ± du rayon issu de la source Si , au 
point de l'espace tridimensionnel M. 
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Pour une onde ultrasonore, 1 •expression des coefficients 
de reflexion Ri est donnee par : 



Ri = 



p 2 c 2 cosy9i -picj 


1 _ C 2 , c 2 „ 2 fl 
. c f cf 


1/2 


P2C2 cos^ + p\c x 


. cf cf 


1/2 



[24] 



10 



15 



10 



:5 



OU : 

- Ci est la vitesse du son dans le milieu M x ; 

- c 2 est la vitesse du son dans le milieu M 2 ; 

- pi est la densite du milieu Mi ; 

- P2 est la density du milieu M 2 . 

Dans cette expression [24] , le terme cosPi pent simplement 
etre estime en fonction des coordonnees dans l'espace du 
point M et du point representant la source Si. 

En se ref grant maintenant a la figure 4B, la m§me 
estimation peut §tre menee pour un point M situe dans le 
milieu M 2 . Dans ce cas f l'onde que recoit le point H est 
une onde transmise par 1' obstacle OBS. En particulier, on 
remarque que les sources de 1- element rayonnant ER ne sont 
plus actives, du fait de 1 • occultation de 1 -element 
rayonnant ER par 1' obstacle OBS. En transmission, le 
raisonnement s> applique comme precedemment avec une 
condition aux limites imposee aux points p a par les 
valeurs des coefficients de transmission T t associes a 
chaque point p t . Dans le cadre de la propagation d'une 
onde ultrasonore, chaque coefficient de transmission T ± 
est donne par la relation : 
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2p 2 Q2 cos fij 



[25] 




1/2 



p 2 o 2 cos Pi +p lCl 1 — I 



Cotnme indique prec€demment , les termes cosPi peuvent etre 



sources S± et du point M. 

Pour estimer les valeurs de sources Si de l 1 obstacle OBS, 
on applique la relation [22] en remplagant toutefois la 
matrice de reflexion R par la matrice de transmission T : 



Dans le cadre de la propagation d'une onde ultrasonore, 
les coefficients des matrices R et T s f estiment pour 
chaque source Si et pour chaque point P A . En particulier, 
chaque coefficient T ifj ou R lfj de la matrice T ou de la 
matrice R (oil i correspond k la i dme ligne et j correspond 
a la j dme colonne) s'exprime en fonction d'un angle Ad 
entre une normale a la surface de l 1 obstacle au point P ± 
et une droite passant par le point p ± et par une source 
Sj. On peut ecrire ainsi, de mani^re g&nerale, les deux 
relations exprimant les valeurs des coefficients des 
matrices T et R par les relations respectives suivantes : 

T U = ft (cos/?ij) [27] 



determines en fonction des coordonnees respectives des 



v = [F(P)]- 1 T[P(P)]v' 



[26] 
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Rij = f r (cos^) [28] 

« 

oil f t est donne par la relation [25] et f r est donne par 
la relation [24] . 

5 

De fagon plus generale, en reference aux figures 4A et 4B, 
on indique que, si l'on considdre 1" obstacle comme un 
materiau plein representant un milieu M2 distinct d'un 
milieu Ml dans lequel se propageait initialement l'onde 

10 principale : 

- pour une reflexion de l'onde principale sur l 1 obstacle 
en tant que milieu M2 (la surface de I 1 obstacle formant 
dioptre entre les milieux Ml et M2) , les hemispheres HEMi 
sont orientes vers l'exterieur de 1' obstacle (figure 4A) ; 

15 - pour une transmission de l'onde principale dans 
1' obstacle, les hemispheres HEMi sont orientes vers 
l'interieur de l 1 obstacle (figure 4B) . 

On se r6f£re maintenant S la figure 5A pour deer ire le cas 
20 d'un l ! obstacle plan OBS de dimensions finies, excite par 
un element rayonnant ER, incline d ! un angle predetermine 
par rapport & I 1 obstacle OBS. Comme indique pr6c€demment , 
pour une onde ultrasonore, on prendra en compte 
1 1 inclinaison de 1 1 element rayonnant pour calculer la 
25 contribution de l'onde emise par l 1 element rayonnant au 
point M. Par ailleurs, de facjon particulierement 
avantageuse, on maille une surface qui englobe la surface 
de l 1 obstacle (figure 5A) . Pour une tranche de I'espace 
delimitee par l 1 element rayonnant, d ! une part, et 
30 l'obstacle, d' autre part (figure 5A) , on peut considerer 
trois types de sources : 
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- les sources S'i de 1 1 element rayonnant ER, 

- des sources SOi, qui renvoient 1 1 onde secondaire, par 
reflexion de l! obstacle OBS, en fonction d'un certain 
coefficient de reflexion R de 1' obstacle ; et 

.5 - des sources SS±, qui ne renvoient pas d'onde secondaire 
et auxquelles un coefficient de reflexion nul peut §tre 
attribu§ si 1' obstacle sepaire deux milieux d' indices 
identiques. Dans ce cas r ces sources SSi sont 
considSrges comrae "Steintes" dans la tranche de 

10 l'espace precitge et ne sont pas prises en compte dans 

les calculs de la grandeur physique au point M de la 
figure 5A. En revanche, ces sources SS± peuvent etre 
actives par reflexion de l'onde principale si 
l 1 obstacle OBS s€pare deux milieux d ! indices 

15 diffgrents. 

Par ailleurs, pour estimer les grandeurs scalaires ou 
vectorielles associees a un point M d'un demi-espace 
delimits par la surface englobant 1 1 obstacle OBS (k droite 
20 de la figure 5B) , on considSre : 

- les sources SO'i de l 1 obstacle, actives par transmission 
de l f onde principale, et 

- les sources SSi, auxquelles on affecte maintenant un 
coefficient de transmission §gal k 1 si 1* obstacle 

25 s§pare deux milieux de memes indices. Ces sources SSi se 

comportent finalement (aux angles d' incidence pres) 
comme les sources S'i de 1' element rayonnant ER. 
Les sources S'i de I'SlSment rayonnant peuvent alors etre 
"eteintes" pour le calcul des grandeurs physiques dans ce 
30 demi-espace. 
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Pour calculer les valeurs v 1 des sources S ! i de l f element 
rayonnant (desquelles sont deduites les valeurs v des 
sources de !■ obstacle selon les xelations [22] et [26] ) , 
on appliquera simplement des conditions aux limites aux 
5 points de la surface active de 1' element rayonnant ER. Par 
exemple, pour une propagation d'onde ultrasonore, on peut 
indiquer que les vitesses sonores aux points de la surface 
de l 1 Element rayonnant ER sont perpendiculaires a cette 
surface et de modules v 0 egaux entre eux. 

10 

De maniere g^ngrale, on indique que l'espace 
tridimensional peut §tre ainsi d^coupe par des interfaces 
deiimitant des milieux de propriety distinctes, chaque 
interface representant un obstacle au sens de la presente 

15 invention. On calcule alors les grandeurs physiques dans 
chaque tranche de l'espace. Par exemple, dans le cadre de 
1 'etude d'une het#rostructure (a plusieurs interfaces), on 
peut appliquer le procede ci-avant pour des tranches 
successives de l'espace en considSrant deux interfaces : 

20 l ? une representant un "Element rayonnant" au sens des 
figures 4A et 5A, par exemple par transmission d'une onde 
regue, et 1 1 autre representant un obstacle recevant l'onde 
transmise. Avantageusement , on tient compte, pour chaque 
tranche de l'espace, des contributions de toutes les 

25 interfaces comme exprime par les relations [10] et [22] . 

Toutefois, dans une realisation pratique preferee, en 
particulier pour programmer la simulation d'une 
interaction, on considerera avantageusement tous les 
30 obstacles de tout l'espace autour d'un point M et l'on 
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imposera une condition sur la position du point M par 
rapport a chaque source pr^sente dans I'espace. 

Pref 6rentiellement , en reference la figure 7A, on teste 

le produit scalaire SMr & chaque iteration par rapport & 
une source S, par exemple sous la forme : 



= cos0 



oH r est le vecteur reliant la source S au point de 
contact P de la demie sphere avec 1 1 element de surface dS 
10 consid6r6, dans le cas oCl l'on pr^voit une seule source 
par hemisphere. Dans le cas oH l'on prevoit plut6t un 

triplet de sources SI, S2 , S3, la base du vecteur r est 
pr€f €rentiellement situSe au barycentre des trois sources 
SI, S2, S3. D'ailleurs, dans le cas de trois sources par 
15 echantillon de surface, le calcul du produit scalaire 
concerne chaque source Si du triplet SI, S2, S3, 

Typiquement, si cos0 est positif, on tient compte de la 
contribution de cette source S dans 1' estimation de 
20 1 1 interaction . 

En revanche, si cos0 est nfigatif, on affecte une valeur 
(scalaire ou vectorielle) nulle & cette source S dans 
1 1 estimation de 1 1 interaction . 



En variante du calcul du produit scalaire ci-avant, on 
peut calculer une "altitude" du point M. Dans ce cas, le 
test porte sur une grandeur du type : 



WO 2004/044790 



PCT/FR2003/003323 



SM.r 




Bien entendu, d'autres types de tests sont possibles. Par 
5 exemple, dans le cas d'un calcul par rapport a 1' angle 0 , 
on peut choisir cet angle dans un c6ne d'ouverture 
choisie, ou autre. 

Finalement, cette demarche permet avantageusement de 
10 systematiser une configuration quelconque des sources par 
rapport au point d 1 observation M, en introduisant 
s implement une §tape supplement aire de test, a chaque 
iteration sur une source S, de la position de cette source 
S par rapport au point M, comme indiqu£ ci-avant. 

15 

Cette demarche s 1 avere particulierement avantageuse pour 
des surfaces k mailler qui sont relativement complexes, 
notamment lorsque le point d' observation M est susceptible 
de se situer dans une zone d ! ombre par rapport & certaines 

20 sources, comme le montre la figure 7B. Sur cette figure 
7B, on a represent^ en traits pointings les demi-sph£res 
associ£es aux sources dans la zone d 1 ombre du point 
d'observation M, et pour lesquelles, par consequent, la 
contribution est f ix£e comme etant nulle dans 

25 1 1 interaction estim€e. Bien entendu, dans le cas ou une 
source est €crant£e par une surface £chantillonn£e, bien 

que le produit scalaire SMr associe & cette source reste 

positif, un second test determine si le vecteur SM 
traverse ou non une surface echantillonn£e . Dans 
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1 'affirmative, cette source est considerge comme inactive 
specif iguement pour la region du point M. 

Ainsi, en termes plus generaux, le precede au sens de 
5 1 ' invention pr^voit pref erentiellement au moins une etape 
supplementaire, pour chaque echantillon de surface, de 
test de la valeur d f un produit scalaire entre : 

- un premier vecteur r normal & 1 1 Echantillon de surface 
et dirige vers le sommet P de 1 'hemi sphere, tel que 

10 represents sur la figure 7A, et 

- un second vecteur SM tire entre une source S associSe a 
cet hemisphere et le point M qui se situe dans la region 
d 1 observation, 

en distinguant, en particulier : 
15 - le cas oil ce produit scalaire est inferieur a un 

seuil predetermine et la contribution de cette source 
n'est pas prise en compte, et 
- le cas ou ce produit scalaire est supErieur £ un 
seuil predetermine et la contribution de cette source 
20 est ef f ectivement prise en compte. 

Dans l'exemple ci-avant oil I'on considere l 1 angle 0 entre 
ces deux vecteurs, le seuil predetermine precite est bien 
entendu la valeur 0 et l'on distingue simplement les cas 
25 ou le produit scalaire est positif ou negatif . 

Bien entendu, ce choix n'est pas limitatif de sorte que 
pour une heterostructure & plusieurs dioptres paralldles, 
on pourra encore considerer, de fagon avantageuse, des 
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semi-espaces successifs, comme d€crit ci-avant en 
reference aux figures 4A et 4B. 

La simulation de la figure 5C correspond, pour une onde 
5 ultrasonore, a la situation des figures 5A et 5B en tenant 
compte : 

- de la contribution de 1- emission de l'onde principale 
par 1' Element rayonnant ER ; 

- de la contribution de la reflexion de cette onde 
10 principale par 1' obstacle ; et 

- de la contribution de la transmission de l'onde 
principale par 1' obstacle. 

Les lignes de niveau de la figure 5C correspondent a 
15 differents paliers de pression acoustique. L' element 
rayonnant ER est dispose a 10 mm de 1> obstacle OBS et 
incline de 20° par rapport a ce dernier. On remarque en 
particulier des f ranges d' interferences dans une zone 
entre 1 » obstacle OBS et 1' element rayonnant ER. Une telle 
20 simulation peut avantageusement indiquer une position 
ideale d'un capteur ultrasonore. Ces capteurs ultrasonores 
comportent habituellement un transducteur comme element 
rayonnant actif et un detecteur pour mesurer les ondes 
ultrasonores recues. La simulation de la figure 5C peut 
25 ainsi indiquer en outre la forme id§ale d'un capteur 
ultrasonore, selon les applications souhaitees, pour une 
forme d' obstacle donn^e. 

La simulation de la figure 5C a ete ef fectuee grace a- un 
30 calcul matriciel programme a l'aide du logiciel de calcul 
MATLAB© . Le nombre de mailles total choisi pour 
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1 ! obstacle et pour 1 1 616ment rayonnant ( ici , quelques 
centaines en tout) est alors optimist : 

d'une part, pour limiter la durfee des calculs ; et 

- d' autre part, pour que la taille des mailles reste 
5 infer ieure k une demie longueur d'onde, de maniere a 

verifier le critdre de Rayleigh. 

On indique toutefois que, comme les elements §l mailler 
dans la mise en oeuvre du proc§d§ selon 1 1 invention sont 
10 simplement des surfaces, les temps de calcul sont bien 
moins longs que ceux necessaires dans la mise en oeuvre 
d'un procede de calcul de type "par elements finis". 

La pr§sente invention peut ainsi se manif ester par la mise 
15 en oeuvre d ! une succession d' instructions d'un produit 
programme d'ordinateur stocke dans la m^moire d ! un disque 
dur ou sur un support amovible et se deroulant comme 
suit : 

- choix d'un pas de maillage notamment en fonction de la 
20 longueur d'onde de I'onde principale ; 

determination des coordonnees des sources Si et/ou S'i 
et des points de contact Pi et/ou P'i ; 

choix d'un type d'onde en jeu et calcul des 
coefficients des matrices d' interaction appliqu^es aux 
25 points Pi et/ou P'i par la mise en oeuvre d'un logiciel 

de calcul matriciel ; 

- choix d'une loi de reflexion et/ou de transmission de 
la surface de 1' obstacle et calcul des coefficients des 
matrices de reflexion et/ou de transmission ; 

30 - calcul des valeurs des sources S± et/ou S'i ; et 
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calcul des grandeurs scalaires ou vectorielles en tout 
point de l'espace tridimensional . 

A ce titre # la pr^sente invention vise aussi un tel 
5 produit programme d 1 ordinateur , stocke dans une mSmoire 
d'unit6 centrale ou sur un support amovible propre a 
coop§rer avec un lecteur de cette unit 6 centrale, et 
comportant en particulier des instructions pour mettre en 
oeuvre le proc€d€ selon l 1 invention. 

10 

Bien entendu, la pr€sente invention ne se limite pas a la 
forme de realisation decrite ci-avant a titre d'exemple ; 
elle s'Stend & d'autres variantes . 

15 Ainsi, on comprendra que, m§me si dans les figures 
commentees ci-avant on repr§sente k la fois la surface 
d'un obstacle et la surface d'un €l§ment rayonnant, la 
pr^sente invention s 1 applique aussi i. 1 1 estimation de 
grandeurs physiques dans le cadre d« une onde interagissant 

20 avec un obstacle et emise en champ lointain. Dans ce 
contexte, il n'est pas necessaire de mat£rialiser la 
surface d'un 616ment rayonnant & mailler et les relations 
[8] et [14] ci-avant suffisent a determiner 1 ' interaction 
entre cette onde et 1 1 obstacle . 

25 

On a indiqu# ci-avant des equations permettant de calculer 
les grandeurs scalaires ou vectorielles en un point M de 
l'espace, pour des ondes electromagnet iques, ou encore 
acoustiques- Bien entendu, ces grandeurs peuvent €tre 
30 estimges pour d ! autres types d 1 ondes, notamment pour des 
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ondes thermiques, des ondes eiectromagnetiques impliquant 
des antennes radiof requences , ou autres . 

Bien entendu, la prSsente invention ne se limite pas & une 
5 application au contrdle non destructif, mais & tout type 
d 1 application, notamment en imagerie medicale, par exemple 
pour 1 ! etude de mi crosyst ernes mettant en oeuvre une 
microscopie acoustique a tniroirs mobiles . 

10 On a decrit ci-avant des interactions entre une onde et un 
obstacle unique. Bien entendu, la presente invention 
s' applique & une interaction avec plusieurs obstacles. A 
cet effet, il convient simplement de mailler les surfaces 
de ces obstacles et d 1 additionner leur contribution pour 

15 1* estimation d'une grandeur vectorielle ou scalaire en un 
point quelconque de l'espace. De mSme, corame indique ci- 
avant, la surface de 1" obstacle OBS peut £tre plane, ou 
encore incurv^e, ou encore de forme complexe quelconque. 

20 Ainsi, dans le cadre d'une onde interagissant avec 
plusieurs obstacles dans l'espace, vine simulation 
Equivalent e a celle representee sur la figure 5C 
permettrait de positionner des capteurs et/ou des Elements 
rayonnants en fonction de la configuration de ces 

25 obstacles, notamment pour une application & la 
determination de la position de hauts-parleurs dans un 
habitacle cloisonne, tel qu'un habitacle de vehicule 
automobile. 

30 L'espace tridimensionnel peut §tre decoupe en une 
pluralite de regions, comme decrit ci-avant en reference 
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aux figures 4A, 4B # 5A et 5B. Toutefois, pour qu'une 
surface de l'une desdites regions soit consid§r6e comme un 
obstacle d'une onde prineipale, devenant actif par 
Emission d'une onde secondaire, 1 1 incidence de l'onde 
5 principale sur cette surface doit pr6f erentiellement 
rester inf 6rieure ou 6gale & 90° . 
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Re vendi c a t i ons 

1. Procede pour evaluer une grandeur physique associ6e a 
une interaction entre une onde et un obstacle, dans une 
5 region de 1 1 espace tridimensional , dans lequel : 

a) on determine, par maillage, une plurality 
d f §chanti lions de surface (dS±) dont une partie au moins 
represente la surface d'un obstacle recevant une onde 
principale et 6mettant, en reponse, une onde secondaire, 

10 et I 1 on attribue & chague £chantillon de surface au moins 
une source (Si) €mettant une onde 61ementaire representant 
une contribution a ladite onde secondaire, 

b) on forme un systeme matriciel comportant : 

- une matrice d' interaction (F (M) ) , inversible, 
15 appliqu^e & une r6gion donn^e (M) de 1 ' espace et 

comportant un nombre de colonnes correspondant & un 
nombre total de sources, 

une premiere matrice colonne dont chaque coefficient 
(Vi) est associg a une source (Si) et caracterise 

20 l'onde glementaire qu'elle 6met, 

et une seconde matrice colonne, obtenue par une 
multiplication de la premiere matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d'une grandeur physique (V(M)) 

25 representative de l'onde Smise par 1' ensemble des 

sources en ladite region donnee (M) , 

c) pour estimer les coefficients de la premiere matrice. 
colonne (v^) , on aff ecte des valeurs de grandeur physique 
choisies (V(Pi)) k des points predetermines (Pi) , associes 

30 chacun & un €chantillon de surface (dSi) , lesdites valeurs 
choisies (V(Pi)) §tant rangges dans la seconde matrice 
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colonne, et I 1 on multiplie cette seconde matrice colonne 
par l 1 inverse de la matrice d 1 interaction appliqu^e 
auxdits points predetermines (Fi) , 

d) pour evaluer ladite grandeur physique (V(M)) 
5 representant l'onde emise par 1' ensemble des sources dans 
une region donnee (M) de I'espace tridimensionnel , on 
applique la matrice d' interaction a ladite region donnee 
(M) et on multiplie cette matrice d' interaction par la 
premiere matrice colonne comportant les coefficients 
10 estimes a l'etape c) . 

2. Procede selon la revendication 1, dans lequel, pour 
^valuer une grandeur physique representative d'une 
interaction entre un Pigment rayonnant une onde principale 
15 et un obstacle recevant cette onde principale, 

- a l'etape a), on determine en outre, par maillage, une 
plurality d'gchantillons de surface (dS'i) representant 
ensemble une surface active de 1" Element rayonnant 1 1 onde 
principale, et l'on attribue & chaque echantillon de la 

20 surface active au moins une source (S ! i) emettant vine onde 
eiementaire representant une contribution I ladite onde 
principale, 

- on applique en outre les etapes b) , c) et d) aux 
echantillons de la surface active, et 

25 - on evalue ladite grandeur physique (V(M)) representant 
1 1 interaction entre l 1 element rayonnant et l 1 obstacle dans 
une region donnee (M) de I'espace tridimensionnel, en 
prenant en compte la contribution, en ladite region donnee 
(M) , de 1 1 onde principale emise par 1 1 ensemble des sources 

30 de la surface active et la contribution de l'onde 



WO 2004/044790 PCT/FR2003/003323 

49 

secondaire 6mise par 1' ensemble des sources de la surface 
de I 1 obstacle. 

3. ProcSde selon l'une des revendi cat ions 1 et 2, dans 
5 lequel chaque coefficient de la matrice d' interaction, 

appliquee §l une region donn6e de l'espace, est 
repr^sentatif d'une interaction entre une source et ladite 
region donn6e et la valeur de chaque coefficient est 
fonction d'une distance entre une source et ladite region 
10 donnee. 

4. Proc£d§ selon l'une des revendications 1 a 3, dans 
lequel la matrice d 1 interaction appliquee, & l f 6tape c) , 
auxdits points pr6d6termin£s (P±) , comporte un nombre de 

15 lignes correspondant a un nombre total de points 
pr6d§termin6s (Pi) . 

5. Proc£d€ selon l'une des revendications la 4, dans 
lequel la grandeur physique ^valuer est une grandeur 

20 scalaire (V(Pj.)) et, a l'€tape a), on attribue une source 
unique a chaque Echantillon de surface. 

6. Proc€d6 selon la revendication 5, dans lequel la 
matrice d' interaction (F (M) ) appliquee, a l»6tape d) , a 

25 une region de 1 1 espace (M) , comporte line ligne. 

7. ProcedS selon l'une des revendications 5 et 6, dans 
lequel chaque point pr6d6termin€ (Pi) associ€ k un 
echantillon de surface (dSi) correspond & un point de 

30 contact entre cet Echantillon de surface (dSi) et -un 
h&nisphdre dont la surface est 6gale & la surface de cet 
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echantillon de surface , et de centre correspondant a une 
position de la source (Si) qui est attribute a cet 
echantillon de surface. 

5 8. Proc6d4 selon l'une des revendications 5^7, dans 
lequel : 

- l'onde principale est une onde glectrique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge glectrique associees chacune a une 

10 source, et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel Slectrique. 

9. Proc£d£ selon l'une des revendications 5 & 7, dans 
15 lequel : 

- l'onde principale est une onde magnet ique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de flux magn£tique associees chacune k une source, 
et 

20 - les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel magnetique. 

10. Proced§ selon l'une des revendications 5 & 7, dans 
lequel ; 

25 - l'onde principale est une onde sonore, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de vitesse de son associees chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
30 valeurs de pression acoustique. 
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11. Proc§d6 selon l'une des revendi cat ions 1 a 4, dans 
lequel la grandeur physique a ^valuer est une grandeur 

;" vectorielle (V(Pi)) exprim€e par ses trois coordonntes 
dans I'espace tridimensional, et l'on attribue, a l'6tape 
5 a), trois sources (SAi, SB if SQl) k chaque 6chantillon de 
surface (dSi) . 

12. Precede selon la revendication 11, dans lequel la 
matrice d' interaction (Fv(M)) appliqu6e, a l'€tape d) , a 

10 une region de I'espace (M) , comporte une ligne pour chaque 
coordonn^e (X,Y,Z) de I'espace. 

13. Proc6d6 selon l'une des revendications 11 et 12, dans 
lequel : 

15. - les trois sources attributes a chaque gchantillon de 
surface sont sensiblement dans un m^me plan et 
- chaque point pr§d§termin6 (Pi) associS & un tchantillon 
de surface (dSi) correspond I un point de contact entre 
cet §chantillon et un hemisphere dont la surface est €gale 

20 k la surface de cet 6chantillon, et de centre 
correspondant & la position d'un barycentre des trois 
sources . 

14. Proctde selon la revendication 13, dans lequel les 
25 trois sources d'un m§me €chantillon de surface forment 
sensiblement un triangle 6quilat6ral, et les triangles des 
gchantillons de surface sont orientes sensiblement 
al^atoirement les uns par rapport aux autres • 

30 15. Proc€d§ selon l'une des revendications 11 3. 14, dans 
lequel : 



WO 2004/044790 



PCT/FR2003/003323 



52 

- l ! onde principale est une onde eiectrique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge eiectrique associees chacune a une 
source, et 

5 - les coefficients de la seconde matrice colonne sont des 
valeurs de coordonn^es d'un champ eiectrique. 

16. Proced§ selon l'une des revendications 11 a 14, dans 
lequel : 

10 - l'onde principale est une onde magnet ique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de flux magnetique associees chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
15 valeurs de coordonnees d ! un champ magnetique. 

17. Procede selon l'une des revendications 11 A 14, dans 
lequel : 

- l'onde principale est une onde sonore, 

20 - les coefficients de la premiere matrice colonne sont des 
valeurs de vitesse de son associees chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de coordonnees d'une vitesse acoustique. 

25 

18. Procede selon l'une des revendications pr€cedentes, 
dans lequel, pour estimer la contribution de I 1 onde 
secondaire en ladite region donnee & 1 1 etape d) , lesdites 
valeurs de grandeur physique (V(Pi)) choisies a 1" etape c) 

30 sont fonction d'un coefficient predetermine de reflexion 
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et/ou de transmission de I'onde principals par chaque 
echantillon de surface de l 1 obstacle, 

19. Procede selon la revendication 18, prise en 
5 combinaison avec l'une des revendications 6 et 12, dans 
lequel l'onde secondaire correspond & une reflexion de 
1 1 onde principale sur 1 1 obstacle et 1 1 hemisphere est 
oriente vers l'exterieur de 1' obstacle. 

10 20. Procede selon la revendication 18, prise en 
combinaison avec l'une des revendications 6 et 12, dans 
lequel l'onde secondaire correspond a une transmission de 
l'onde principale dans 1' obstacle et 1' hemisphere est 
oriente vers 1 ' interieur de l f obstacle. 

21. Procede selon l'une des revendications 19 et 20, dans 
lequel, k 1 1 etape c) , on determine les valeurs (v'i) 
associees aux sources (S'i) de l 1 element rayonnant (ER) et 
on met en forme au mo ins : 
20 - une premiere matrice d' interaction (F(P)) representant 
la contribution des sources de l 1 obstacle aux points 
predetermines de la surface de 1' obstacle (Pi), 

- une seconde matrice d' interaction (P f (P)) representant 
la contribution des sources de l 1 element rayonnant aux 

25 points predetermines de la surface de l 1 obstacle (Pi), 

- une matrice de reflexion (R) ou de transmission (T) , 
dont les coefficients representent des coefficients de 
reflexion ou de transmission en chaque point predetermine 
(Pi) de 1 'obstacle, 

30 pour determiner les valeurs des sources de 1' obstacle (vi) 
en f onction des valeurs des sources de 1 1 element rayonnant 
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(v'i) et d'une multiplication des premiere et seconde 
matrices d» interaction et de la matrice de reflexion ou de 
transmission. 

22. Procedg selon la revendication 21, dans lequel, a 
l«6tape c) , on determine les valeurs (v'i) associSes aux 
sources (S'i) de l 1 Element rayonnant (ER) en tenant compte 
de la reception de l'onde secondaire par 1' element 
rayonnant (ER) et en mettant en forme en outre : 
- tine troisiSme matrice d' interaction (F(P')) reprSsentant 
la contribution des sources de l 1 obstacle aux points 
predetermines de la surface de 1» element rayonnant (P'i) , 

et une quatridme matrice d' interaction (F» (P 1 ) ) 
repr€sentant la contribution des sources de 1 1 element 
rayonnant aux points predetermines de la surface de 
1 1 e 1 ement rayonnant ( P 1 ± ) . 

23. Proc€d6 selon l'une des revendi cat ions 19 a 22 , dans 
lequel la surface de 1" obstacle correspond a une interface 
entre deux milieux distincts d'une hetSrostructure . 

24. Procede selon l'une des revendi cat ions pr6c6dentes, 
dans lequel l'onde principale est une onde sonore et les 
coefficients de la matrice d' interaction sont chacun 
fonction d'un angle d 1 incidence d'une onde eiementaire 
issue d'une source en ladite region donnee (M) . 

25. Procede selon l'une des revendications 7 et 13, dans 
lequel, pour chaque €chantillon de surface, on teste la 

30 valeur d'un produit scalaire entre : 
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- un premier vecteur (r) normal a 1 1 echantillon de surface 
et dirige vers le sommet (P) de 1 'hemisphere (fig. 7A) ,' et 

- un second vecteur ( SM ) tire entre une source (S) 
associee a cet hemisphere et ladite region donn£e (M) , 

5 en distinguant : 

- le cas ou ce produit scalaire est infer ieur a un 
seuil predetermine et la contribution de cette source 
n'est pas prise en compte, et 

- le cas oil ce produit scalaire est sup£rieur a un 
10 seuil predetermine et la contribution de cette source 

est ef f ectivement prise en compte. 

26. Procede selon l'une des revendications precedentes, 
dans lequel l'onde principale est une onde sonore et, a 
15 l'etape a)/ on choisit un nombre total d 1 echanti lions de 
surface (dSj.) sensiblement en fonction d'une longueur 
d'onde de I'onde sonore pour verifier le critdre, de 
Rayleigh. 

20 27. Procede selon l'une des revendications precedentes, 
dans lequel on compare une pluralite de valeurs de la 
grandeur physique estimee k 1 1 etape d) , obtenues pour une 
pluralite de regions de l'espace, pour select ionner une 
region candidate pour la disposition d'un element 

25 rayonnant destine k interagir avec l 1 obstacle. 

28. Procede selon l'une des revendications 2 a 27, dans 
lequel 1 1 element rayonnant est un capteur, pour un 
contr61e non destructif , destine & analyser un objet 
30 formant un obstacle de l'onde principale. 
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\ 

29. Produit programme d'ordinateur, stocks dans une 
m^moire d ! unite centrale ou sur un support amovible propre 
a coop^rer avec urs lecteur do cette units centrale, 
caracterise en ce qu'il comporte des instructions pour 
5 mettre en ceuvre le procedS selon l'une des revendications 
precedentes . 



WO 2004/044790 



PCT/FR2003/003323 




WO 2004/044790 



PCTYFR2003/003323 



2/9 




WO 2004/044790 



PCT/FR2003/003323 



3/9 



SB2 




FIG. 2C 



WO 2004/044790 



PCT/FR2003/003323 




FIG. 4B 



WO 2004/044790 



PCT/FR2003/003323 



6/9 




FIG. 5B 



WO 2004/044790 



PCT/FR2003/003323 



7/9 




FIG. 5C 



WO 2004/044790 



PCT/FR2003/003323 



8/9 



FIG. 6 




WO 2004/044790 



PCT/FR2003/003323 




FIG. 7B 



INTERNATIONAL SEARCH REPORT 



Intematit plication No 

PCTAk ui/03323 



A. CLASSIFICATION OF SUBJECT MATTER 

PC 7 G06F17/50 



According to international Patent Classification ((PC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 G06F 



Documentation searched other than minimum documentation to the extent thai such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 

EPO-Internal , INSPEC, COMPENDEX 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category • CilaUon of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



LIEBEAUX N ET AL: "The distributed 

sources method: a concept for open 

magnetic cores modelling" 

EUROPEAN PHYSICAL JOURNAL, APPLIED 

PHYSICS, NOV. 2002, EDP SCIENCES, FRANCE, 

'Online! 

vol. 20, no. 2, 

25 October 2002 (2002-10-25), pages 

145-150, XP008021970 

ISSN: 1286-0042 

Retrieved from the Internet: 

URL : http : //www . edpsd ences . org/art 1 cl e s/ep 

j ap/pdf /2002/1 l/ap01142 . pdf > 

'retrieved on 2003-09-16! 

the whole document 

-/- 



1-10,28, 
29 



11-27 



HI 



Further documents are listed In the continuation of box C. 



□ 



Patent family members are listed in annex. 



° Special categories of cited documents : 

"A' document defining the general state of the art which Is not 

considered to be of particular relevance 
"E" earlier document but published on or after the international 

filing date 

*L' document which may throw doubts on priority claim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

'O' document referring to an oral disclosure, use, exhibition or 
other means 

"P" document published prior to the international filing date but 
later than the priority date claimed 



"T* later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
Invention 

•X 1 document of particular relevance; the claimed Invention 
cannot be considered novel or cannot be considered to 
Involve an Inventive step when the document Is taken alone 

"Y* document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person j 
In the art 

document member of the same patent family 



Date of the actual completion of the international search 



15 March 2004 



Date of mailing of the International search report 



25/03/2004 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patent laan 2 
NL-2280HVRl]swijk 
Tel (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax: (+31-70)340-3016 



Authorized officer 



Gulngale, A 



Forni PCT7ISA/210 (second sheet) (July 1992) 



INTERNATIONAL SEARCH REPORT 



Intematio plication No 

PCT/Fiv ui/03323 



C(ConUnuaUon) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ' 



Citation of document, with Indlcatlon.where appropriate, of the relevant passages 



Relevant to claim No. 



PLACKO D ET AL: "Theoretical study of 
magnetic and ultrasonic sensors: 
dependence of magnetic potential and 
acoustic pressure on the sensor geometry" 
ADVANCED NONDESTRUCTIVE EVALUATION FOR 
STRUCTURAL AND BIOLOGICAL HEALTH 
MONITORING, NEWPORT BEACH, CA, USA, 6-8 
MARCH 2001, 

vol. 4335, 6 March 2001 (2001-03-06), 
pages 52-62, XP008021972 
Proceedings of the SPIE - The 
International Society for Optical 
Engineering, 2001, SPIE-Int. Sbc. Opt. 
Eng, USA 
ISSN: 0277-786X 
the whole document 



1-10,28, 
29 



PLACKO D ET AL: "Theoretical computation 
of acoustic pressure generated by 
ultrasonic sensors in the presence of an 
interface" 

SMART NONDESTRUCTIVE EVALUATION FOR HEALTH 
MONITORING OF STRUCTURAL AND BIOLOGICAL 
SYSTEMS, NEWPORT BEACH, USA, 18-20 MARCH 
2002, 

vol. 4702, 18 March 2002 (2002-03-18), 
pages 157-168, XP008021971 
Proceedings of the SPIE - The 
International Society for Optical 
Engineering, 2002, SPIE-Int. Soc. Opt. 
Eng, USA 
ISSN: 0277-786X 
the whole document 



11-27 
1-29 



Form PCT7ISA/210 (continuation ol second sheet) (July 1992) 



RAPPORT DE RECHERCHE INTERNATIONALE 



Oemande donate No 

PCT/FR'D3/03323 



A. CLASSEMENT DE L'OBJET DE LA DEMANDE 

CIB 7 G06F 17/50 



Selon la classification Internationale des brevets (CIB) ou a la fols seton la classification nationals at la CIS 

B. DOMAINES SUR LESQUELS LA RECHERCHE A PORTE 

Documentation mlnlmalo consultee (systems de classification suM des symboles de dassemenl) 

CIB 7 S06F 



Documentation consultee autre que la documentation mlnimale dans la mesure oil ces documents relevenl des domahes sur lesquels a ports la rechsrchs 
Bass de donnees slectronique consultee au cours de la recherchs IntsrnaUonals (nom de la bass ds donnees, st si realisable termes ds recherche utilises) 

EPO-Internal, INSPEC, COMPENDEX 



C. DOCUMENTS CONSIDERES COMME PERTINENTS 



Categorte ° Idsntlflcatlon dss documents cites, avec. Is cas echeant. rindlcatton des passages pertinents 



no. des revendlcations vlsees 



LIEBEAUX N ET AL: "The distributed 
sources method: a concept for open 
magnetic cores modelling" 
EUROPEAN PHYSICAL JOURNAL, APPLIED 
PHYSICS, NOV. 2002, EDP SCIENCES, FRANCE, 
'Online! 

vol. 20, no. 2, 

25 octobre 2002 (2002-10-25), pages 
145-150, XP008021970 
ISSN: 1286-0042 
Extrait de 1 'Internet: 

URL : http : //www . edpsd ences . org/art 1 cl es/ ep 
j ap/pdf/2002/1 l/ap01142 . pdf > 
'extrait le 2003-09-16! 
le document en entler 



1-10,28, 
29 



11-27 



-/-- 



[x] v <* la 



suite du cadre C pour la fin de la llste des documents 



□ 



Lbs documents de families de brevets sont Indlques en annexe 



• Categories speclales de documents cites: 

•A' document deflnlssant Tetat general de la technique, non 

considers comma partlcuUerement pertinent 
*E a document anterisur, mate publle a la date de depot international 

ou apres cette dale 
*L" document pouvant jeter un doute sur une revendication de 
priority ou cite pour determiner la date de publication d'une 
autre citation ou pour une raison specials (telle qulndlquee) 
■O* document se referant a une divulgation orale, a un usage, a 

une exposition ou tous autres moyens 
■p* document publle avant la date de depot international, mais 
posterieurement a la date de priortte revendiquee 



*V document ulterieur publle apres la date de depot International ou la 
date de priorite et n'appartenenant pas a fetat de la 
technique pertinent, mais cite pour comprendre le principe 
ou latheorie constituent la base de Hnventlon 

•X' document particullerement pertinent; rinven tlon revendiquee ne peut 
etre consideree comme nouvelle ou comme impllquant une acttvtte 
Inventive par rapport au document considers Isolement 

•V document particullerement pertinent; rinven tlon revendiquee 

ne peut etre consideree comme impllquant une activtt6 Inventive 
lorsque le document est associe a un ou plusieurs autres 
documents de meme nature, cette comblnaison etant evldente 
pour une personne du metier 

document qui fait partle de la meme famine de brevets 



Date a laqueUe la recherche Internationale a ete effectivement achevee 



15 mars 2004 



Date d'expeditlon du present rapport de recherche Internationale 

25/03/2004 



I Nom et adresse postale de radmlnistration charges de la recherche Internationale 
Office Europeen des Brevets, P.B. 5818 Patentlaan 2 
NL-2280HVRqswiik 
Tel (+31-70) 340-2040. Tx. 31 651 epo nl, 
Fax: (431-70) 340-3016 



Fonctlonnaire autorise 



Gulngale, A 



«^.(tLo\ /1,»UU4 tQQO\ 



RAPPORT DE RECHERCHE INTERNATIONALE 



Demands lj rationale No 

PCT/FR .3/03323 



C(sulte) DOCUMENTS CONS1DERES COMME PERTINENTS 



Cat£gorfe ° Identification des documents cites, avec.le cas echeant, I'indlcationdes passages pertinents 



no. des revendicaUons visees 



PLACKO D ET AL: "Theoretical study of 
magnetic and ultrasonic sensors: 
dependence of magnetic potential and 
acoustic pressure on the sensor geometry" 
ADVANCED NONDESTRUCTIVE EVALUATION FOR 
STRUCTURAL AND BIOLOGICAL HEALTH 
MONITORING, NEWPORT BEACH, CA, USA, 6-8 
MARCH 2001, 

vol. 4335, 6 mars 2001 (2001-03-06), pages 

52-62, XP008021972 

Proceedings of the SPIE - The 

International Society for Optical 

Engineering, 2001, SPIE-Int. Soc. Opt. 

Eng, USA 

ISSN: 0277-786X 

le document en entier 

PLACKO D ET AL: "Theoretical computation 
of acoustic pressure generated by 
ultrasonic sensors 1n the presence of an 
interface" 

SMART NONDESTRUCTIVE EVALUATION FOR HEALTH 
MONITORING OF STRUCTURAL AND BIOLOGICAL 
SYSTEMS, NEWPORT BEACH, USA, 18-20 MARCH 
2002, 

vol. 4702, 18 mars 2002 (2002-03-18), 

pages 157-168, XP008021971 

Proceedings of the SPIE - The 

International Society for Optical 

Engineering, 2002, SPIE-Int. Soc. Opt. 

Eng, USA 

ISSN: 0277-786X 

le document en entier 



1-10,28, 
29 



11-27 
1-29 



Foimubiro PCT/1SA/21 0 (suite do la deuxfeme feuMa) (juilet 1892) 



